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Abstract. While considering a class of generalized negative binomial 
£h states, we verify that the basic minimum properties for these states to be 

considered as coherent states are satisfied. We particularize them for the 
case of the Hamiltonian of the isotonic oscillator and we determine the 
corresponding Husimi's Q-function. This function may be used to deter- 
mine a lower bound for the thermodynamical potential of the Hamilton- 
ian by applying a Berezin-Lieb inequality. 

Ph 1 Introduction 

^ The negative binomial states (NBS) are the field states that are super- 

position of the number states with appropriately chosen coefficients JH. 
Precisely, these labeling coefficients are such that the associated photon- 
counting distribution is a negative binomial probability distribution |0. 

t— I As matter of fact, these coefficients turn out to be an orthonormal basis of 

a weighted Bergman space of analytic functions on the complex unit disk 
satisfying a certain growth condition. Furthermore, the NBS are considered 
as intermediate states between pure coherent states and pure thermal states 

^ and reduce to Susskind-Glogower phases states for a particular limit of 

the parameter 01 . 

Now, as in J5j |6|, we replace the labeling coefficients by an orthonormal 
basis of a Hilbert space that generalize the weighted Bergman space of an- 
alytic functions on the unit disk we have mentioned to consider a class of 
. — generalized negative binomial states (GNBSs) in this sense. Here, we pre- 

^ cisely verify that the basic minimum properties for the constructed states 

to be considered as coherent states are satisfied. Namely, the conditions 
which have been formulated by Klauder 0: (a) the continuity of labeling, 
(b) the fact that these states are normalizable but not orthogonal and (c) 
these states fulfilled the resolution of the identity with a positive weight 
function. 

Next, we particularize the GNBSs formalism for the case of the isotonic 
oscillator (IO) |8| whose importance consists in the fact that it admits ex- 
act analytic solutions an being in a certain sense an intermediate potential 
between the three dimensional harmonic oscillator potential and other an- 
harmonic potentials such as Poschl-Teller or Morse potentials [9J. Next, we 
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determine the Husimi's Q-f unction [10] which turns out to be the expec- 
tation value in the GNBSs representation of the heat semigroup operator 
associated with the IO. The obtained Q-function can be considered as a 
lower symbol for this heat semigroup operator. Finally, a lower bound for 
the thermodynamical potential of the IO may be obtained by applying a 
Berezin-Lieb inequality | TT][T2| . 

The paper is organized as follows. In Section 2, we review briefly the co- 
herent states formalism we will be using. Section 3 deals with some needed 
facts on the generalized weighted Bergman spaces on the disk. In Section 
4, we attach to each of these spaces a set of coherent states generalizing 
the negative binomial states and verify that they satisfy the basic minimum 
properties of coherent states. In Section 5, we recall briefly some needed 
spectral properties of the isotonic oscillator Hamiltonian . In Section 6, we 
obtain the Husimi's Q-function associated with the IO in the coherent state 
representation and we deduce a lower bound for the thermodynamical po- 
tential of the IO. 



2 Coherent states and Berezin-Lieb inequalities 

Here, we review a coherent states formalism starting from a measure space 
"it as a set of data" as presented in [13]. Let X = {x \ x G X} be a set 
equipped with a measure d\i and L 2 (X,d}i) the space of df/-square inte- 
grable functions on X. Let A 2 C L 2 (X,d}i) be a subspace of infinite di- 
mension with an orthonormal basis Let % be another (functional) 

space with dim'H = oo and {(pj}°°_^ is a given orthonormal basis of %. 

Then consider the family of states {| x >} xeX in 7i, through the following 
linear superpositions: 



i 



*>:= (AT (*))-* \<Pj> (2-1) 

/=0 



where 



N{x) = Y J <*>j (*)<*>; (*)• (2-2) 

These coherent states obey the normalization condition 

(x \x) n = l (2.3) 
and the following resolution of the identity of 7i 



I | x >< x | M (x) dfi (x) (2.4) 

which is expressed in terms of Dirac's bra-ket notation | x >< x \ meaning 
the rank-one -operator cp 1-4 (<p \ x)^. | x >. The choice of the Hilbert 
space T-L define in fact a quantization of the space X by the coherent states 



in (2.1 1, via the inclusion map X 3 x h- >|x>G'H and the property (2.4 1 
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is crucial in setting the bridge between the classical and the quantum me- 
chanics. 

Now, given a set of coherent states {| x >}, the concept of upper and 
lower symbols of an operator A was separately introduced by Berezin fTTfl 
and Lieb [12] by 

A= / dji(x)A\x><x\ (2.5) 



to define the upper symbol A of A, and the expectation value 

A (x) := (x | A | x) (2.6) 

for the definition of the lower symbol A of A. Note that given an operator A 
its upper symbol is not unique in general. It can be proved |T2) that given 
any convex function (p, the following inequalities 

j^{A)d ¥ {x)<Tr{ ( p{A))< j^{A)d V {x) (2.7) 

hold and are called Berezin-Lieb inequalities. 

3 Generalized Bergman spaces on D 

Let D = {z G C, |z| < 1} be unit disk endowed with its usual Khaler 
metric ds 2 = —ddLog (1 — zz) dz <8> dz. The Bergman distance on D is given 

,2,/ \ (1 — Zw)(l — Zw) .„„. 

cosh 2 d z, w) = j- 4 3.1 

(1 — zz) (1 — WW) 

and the volume element reads 

1 

with the Lebesgue measure dv (z) . Let us consider the 1— form on E) de- 
fined by 9 = —i (^d — d^j Log (1 — zz) to which the Schrodinger operator 

H,:= (d + i^ext{9)y (d + i^ext {0)) (3.3) 

can be associated. Here a > is a fixed number, d denotes the usual ex- 
terior derivative on differential forms on E) and ext (9) is the exterior mul- 
tiplication by 9 while the symbol * stands for the adjoint operator with 
respect to the Hermitian scalar product induced by the Bergman metric ds 2 
on differential forms. Actually, the operator H a is acting on the Hilbert 
space L 2 (D, dji (z)) and can be unitarly intertwined as 

(l-zz)^A a (l-zz)-^ = H a (3.4) 
in terms of the second order differential operator 

^:=-4(l-zz)((l-zz)Jg-<4). <3 ' 5) 



dy, (z) = —^dv (z) (3.2) 
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The latter one is acting on the Hilbert space 

L 2,cr (D) = L 2 (d, (1 - zzf- 2 dv (z) ) . (3.6) 



The spectral analysis of A a have been studied by many authors, see 
and references therein. Note that this operator is an elliptic densely de- 
fined operator on L 2,a (ID) and admits a unique self-adjoint realization that 
we denote also by A^. The part of its spectrum is not empty if and only 
if a > 1. This discrete part consists of eigenvalues occurring with infinite 
multiplicities and having the expression 

:= Am (a - m - 1) (3.7) 

for varying m = 0, 1, [(a — l)/2] . Here, [x] denotes the greatest integer 
not exceeding x. Moreover, it is well known that the functions given in 
terms of Jacobi polynomials l[T5llT7t by 



cB^M- / (c — 2m — 1) klT (a — m) 

** (Z) "V 7zm\T(a-2m+k) M 

X (-if (1 - ZZ)-" 1 ^n-kp(m-kA-2m-l) ^ _ ^ 

constitute an orthonormal basis of the eigenspace 

Alf (D) := {<p G L 2,a (D) , K<p = e a m( p} . (3.9) 

of A^ associated with the eigenvalue in ( |3.7[ ). Finally, the L 2 — eigenspace 
Al' a (D) = {<p e L 2 ' a (D),A a (p = 0} corresponding to m = and asso- 
ciated to Cq = in ( |3.7[ | reduces further to the weighted Bergman space 
consisting of holomorphic functions <p : D —¥ C with the growth condition 

[ \d> (z) I 2 (1 - zzf- 2 dv (z) < +oo. (3.10) 
Jd 



This is why the eigenspaces in ( |3.9| l have been called generalized Bergman 
spaces on the complex unit disk. 

Remark 3.1. In |5l, we have used the basis in ( |3.8| l to perform a class of coherent 
states belonging to the Hilbert space L 2 (R+, x~ J dx) . The associated coherent 
state transform have been considered as a generalization of the second Bargmann 
transform (HE p.203]). 

4 Generalized negative binomial states 

The negative binomial states are labeled by points z G ID and are of the 
form 




z /£ r,0 >:= (1 - zzY a V W =^ , , V | ^ >, (4.1) 



where a > 1 is a fixed parameter and | ipi > are Fock states. Their photon 
probability distribution 

|< ft | z,a >| 2 = (1 - zzf («)' (4.2) 

obeys the negative binomial probability distribution with parameters A = 
zz and c. Also, observe that the coefficients in the superposition ( |4.1) : 



(43) 

constitute an orthonormal basis of the eigenspace ^Iq' 17 (D) associated with 
the first eigenvalue €q = and consisting of analytic functions on ID with 
the growth condition ( 3. 1Q| > . 

For instance, let a > 1 and m = 0, 1, [{a — 1) /2] be fixed parameters 
and let { | xpk > }j*L be a set of Fock states in a Hilbert space T-L. Then a class 
of generalized negative binomial states (GNBS) can be defined as in OHJ 
by 

\z,<r,m >= [Ma,m (z))~ 2 E °r ( z ) I ft >/ ( 4 - 4 ) 

*c=0 

where A/^m (z) is a normalization factor and {O^'"' (z) is the orthonor- 
mal basis of the generalized Bergman space in \3 .9) . Now, one of the impor- 
tant task to do is to determine is the overlap relation between two GNBSs. 

Proposition 4.1. Let a > 1 and m = 0,1,..., [(c— 1) /2] . T/zen, for every 
z,w € ID, the overlap relation between two GNBSs is given through the scalar 
product 



(a -2m- l)r (a - m) (Af (z) N (w))~ z 
Tzml (—1) T (a — 2m) (1 — zw) 
/ (1 — zw) (1 — wz)\ m ( _ (l — zz)(l — ww) 



(1 — zz) (1 — w) / \ ' ' ' (1 — zw) (1 — aJz) 

w/zere 2F\isa terminating Gauss hypergeometric sum. 

Proof. In view of Eq. ( |4.4| l, the scalar product of two GNBS | z ,c~,m > and 
I w ,o~,m > in H reads 



< w,a,m \z,(r,m > n = (N a , m (z) A/^m (a?)) 2 £ °5t' m ( z ) ( w ) 



(c--2m-l)(Af(z)Af(w))-i 

n{(l-zz){l-ww)) m z ' w ' K ' 
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where 

r^-m)^)"' +^ g / i \ k 

m\T(a-2m) ^ (a - 2m) k \zw ) K ' 

x (x _ 2z z) p^-^-^-D (i _ 2a w) . 

Making use of the following identity ([19. p. 1329]): 

E tt^^™' 00 ^ n ' a) GO = ( i - z - 1) (y - 1) t 



(4.8) 

for n = k,t = 1/zw, j = m, a = a — 2m — 1, x = 1 — 2zz and y = 1 — 2w w, 
we obtain, after calculations, the expression 

«r,m = r(o--m) (-l) m ((l-zw) (l-a;z)) m 
z ' w m!r((7-2m) (l-zw) 17 ^ J 

/ (1 — zz) (1 — ww) \ 

x 2 Fi -m,a-m,a-2m;j- — — ^ . 

\ (1 — zw) (1 — wz) J 

Returning back to Eq. (4.6 1 and inserting the expression ( |4.9| l we arrive at 
the announced formula. ■ 

Corollary 4.2. The normalization factor in ( |4.4| | is given by 

k r / \ (c — 2m — 1) 

AU z = ^ (4.10) 

7T (1 — ZZ) 

for every z G D. 

Proof. We first make appeal to the relation (||T7l p. 212]): 

2 Fi f -n, n + k + Q + 1, 1 + k; ^ = f 1 + , K) , (r) (4.11) 



2 / r(l+K + n 

connecting the 2F 1— sum with the Jacobi polynomial for the parameters 
n = m, K = o~ — 2m — 1, q = and the variable 



T = j _ 2 (1-ZZ)(1-IDW) 

(1 — zw) (1 — wz) 



(4.12) 



to rewrite Eq. ( |4.5| l as 

V^(z)^ m H = (-irv-^-mi-^r (413) 

5T < w; cr, m \ z , cr, m 

x p (cr-2m-l,0) A _ 2 (1 - zg (1 - WW) 

\ (1 — zw) (1 — wz) 
The factor A/"^ m (z) should be such that 

< z,o~,m | z,cr,m >%= 1. (4.14) 



So that we put z = w in ( 4.13[ > and we use the symmetry identity ( |T7l 
p.210]): 

p4 T ' e) (?) = (-ir* 7) (-o (4.i5) 

to obtain the expression 

Ma,m (Z) = {a - 2m Z^ P^- lm - l) (1) (4.16) 
71 (1 — ZZJ 

Finally we apply the fact that (HZl p.1329]): 

M ( ) £(n + a + l) 

in the case of a = 0, n = m and q = a — 2m — 1. This ends the proof. ■ 

Proposition 4.3. Let a > 1 and m = 0,1, [(a — 1) /2] . Then, the GNBS in 
( |4.4| | satisfy the following resolution of the identity 

l n = | z,a,m >< z,a,m \ du um (z) , (4.18) 
Jd 

zffoere 1% is f/ze identity operator and d}i a/m (z) zs a measure which can be ex- 
pressed through a Meijer's G-function as 

o™ 1 Wi ' - . -1 



d M(7 , m (z) := 7T- > - 2m - 1)G{{ I -zz \ ^ J dv (z) , (4.19) 

and dv (z) being- iize Lebesgue measure on D. 

Proof. Let us assume that the measure takes the form 

(z) = M a , m (z) O (z) dv (z) , (4.20) 

where O (z) is an auxiliary density to be determined. Let q> G % and let us 
start by writing the following action 

O [q>] := fj^ \z,a,m >< z,a,m\ d]i UiM (z) \ [cp] (4.21) 

= / < (p | z,a,m >< z,a,m \ d]i am (z) . (4.22) 
Jd 

Making use Eq. §4.4) , we obtain successively 

= / < <p | (A/^m (z)) z X] °fc' m ( z ) I #t >>< z,a,m\ du a , m (z) 

(4.23) 

n +00 ^ 1 

= / E °fc' m ( z ) < ? I >< z ' a ' m I ( z )) 1 d ^> m ( z ) 



o 



(4.24) 



(4.25) 
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We replace the measure dy.ff,m (z) by the expression in the right hand side 
of Eq. ggU) , then Eq. ( |425| can be written without <p as follows 



0= E 

j,k=0 



ID 

Therefore, we need to have 



; ^( Z )0^ (z) Q (z) dv (z) 



>< ip k | . (4.26) 



| d <Ef» (z)O^ (z) Q (z) dv (z) = fy. (4.27) 

For this we recall the orthogonality relation of the 0^' m (z) in the Hilbert 
space L 2 ' a (D), which reads 

JD ' 



This suggests us to set O (z) := (1 — zz) a ~ 2 . Therefore, we get that 

(a -2m- 1) 

71 (1 — ZZ 

By making us of the identity 



d H<r,m (z) = ^ -^dv (z) . (4.29) 



11 /> I a 



J ll 



n : =r(i- fl +fc)c*(i+o ( 4 - 3 °) 



for £ = — zz, a = —1 and b = 0, we arrive at the expression of the measure 
in ( 4.19p . Therefore, Eq. ( 4.26 1 reduces to 

O = E \ti>< Vfc 1= **■ (4-31) 

;,it=0 

The proof is finished. ■ 

Proposition 4.4. Let cr > 1 and m = 0,1,..., [((7 — 1) /2] . Then, the states 
| z, c, m > satisfy the continuity property with respect to the label z G ID. That is, 
the norm of the difference of two states 

d a ,m (z,w) : = ||(| z,a,m > - \ w,a,m >)\\ u (4.32) 

goes to zero whenever z — > w. 

Proof. By using the fact that any GNBS is normalized by the factor given 
in ( 4.10| |, direct calculations enable us to write the square of the quantity in 
14.32) as 

d 2 m (z,w) = 2 (1 — < z, cr,m \ w,a, m >) . (4.33) 

Next, we use of the expression of the scalar product in ( |3.9| l form which 
it is clear that the overlap takes the value 1 as z — > w and consequently 
dcr /m (z,w) — > 0. ■ 

We end this section by the following remarks. 
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Remark 4.5. By a general fact on reproducing kernels [21 J, the proof of proposition 



2.1) also says that the knowledge of the explicit orthonormal basis in ( |3.8[ ) leads 



directly to expression of the reproducing kernel of the generalized Bergman space 



in ( |3.9| via calculations using the formula ( |4.8[ ) due to A. Srivastava and A. B. 
Raolm. 



Remark 4.6. In |6]|, we have used the same basis ( |3.8| under another form as la- 
beling coefficients in order to consider the photon-counting probability distribution 
with the mass function 

ft' m ' X) ■■= 7<r,m,k (1 " Af- 2 "' Al-*l (pf m Z\*-^ (! - 2A )) 2 ' ( 434 ) 

where k = 0,1,2,..., j aiM/ k > is a constant, m = 0,1, [a — 1)/2] and 
A = zz. We have calculated the associated Mandel parameter [22] and we have 
discussed the classicality /nonclassicality of the GNBS with respect to the location 
of their labeling points z inside the hyperbolic disk D. Similar results, in the Eu- 
clidean plane and the Riemann sphere settings have been obtained respectively in 
H and EH. 



Remark 4.7. The fact that we have written the measure d}i aim (z) in ( |4.19 l in 
terms of the Meijer's G-f unction could be of help when tackling the "photon-added 
coherent states (PACS)" problem for the GNBS under consideration. 

5 The isotonic oscillator L a 

Not all quantum Hamiltonians are known to have exact solutions. An 
important model of a solvable class is the isotonic oscillator [25 1 




L « : =d-y^+ ;t +- 7 — - ,«>l/2 (5.1) 



acting in the Hilbert space V. := L 2 (R+, dx) and the eigenf unctions ip G % 
satisfy the Dirichlet boundary condition tp (0) = 0. This operator appears 
in the literature under many names such as Gol'dman-Krivchenkov Hamil- 
tonian [|26| or pseudoharmonic oscillator (8l or Laguerre operator f27\ . It 
is the generalization the harmonic oscillator in three dimensions where the 
generalization lies in the parameter v = a 2 — 1/4 ranging over (0, +oo) 
instead of the angular momentum quantum numbers / = 0,1,2,... . This 
operator may be factorized as follows 

L a = D+D- + oc + 1 (5.2) 

in terms of the operators, having the form 



(5.3) 

It is well known that the Hamiltonian L a admits exact solutions of the form 

r k (x) = ( r^+^+i) ) 1 (* 2 ) . (5-4) 
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where (.) denotes the La guerre polynomial Ifl6l and the corresponding 
eigenvalues are given by ( 112811 ): 

kl = 2a + k + l,k = 0,1,2,... . (5.5) 

Note also that the functions can be obtained by A;— fold application of a 
creation opera tor to the ground state wavef unction ipfi. The vectors {| ip% >}^ 
satisfy the orthogonality relation 

r+oo 

<V£ I r k )n = / # (*) (*) dx = ( 5 - 6 ) 

J o 

and constitute a complete orthonormal basis for the Hilbert space TL. Fur- 



thermore, they can be used together with the eigenvalues in (5.5 1 to define 
the heat semigroup associated with L a as 

e -*«\f);=£e- a i<f\r k >n-r k (5-7) 

for any function / 6 L 2 (R+,dx) .It is also well known that by using the 
Hille-Hardy formula ( fT71 p.242]), this semigroup has an integral represen- 
tation, i.e., 



where 



e- tt *\M*)= W t (x,y)f(y)dy r (5.8) 
Jo 

Here l K (.) denotes the modified Bessel function of the first kind and order 
*(E3p.66]). 

6 HUSIMI'S Q-FUNCTION ATTACHED TO L K 

For a > 1 and m = 0,1, [(c — 1) /2] . A class of generalized negative 
binomial states (GNBS) attached to the isotonic oscillator L a can be defined 
by setting 

|z,<r,m,« >= (A/cr,m (z))" z ( z ) I ^ > 

k=0 



where A/"^ (z) is the factor in JPo) , (z)} are defined in ( |3.8| l and 



{ | t/>£ >} are the Fock vectors given in <\5A) . The diagonal representation 
of e~ tL « in the GNBSs representation is now precised as follows. 

Definition 6.1. The Husimi's Q-function attached to the operator L a is given 
through the mean value 

Q t m (L K )(z) =E { | 2/W>} (e- tK ) =<z,<T,m,u\e- tK \z,a,m,a> (6.2) 

with respect to the set of GNBSs defined in ( |6.1| . 
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Proposition 6.2. The mean value defined in ( |6.2| has the following expression 

_ njl-zzfe-^ ( (zz-e^]^zze^)\ '" 
{I ™ >}1 j " (a-2m-l) { (1 _^) 2 ) 

1 " Z " Y P ( :- lm ' Q) ( 1 + f - 2e ^!~ Z ~ Z) -^ ) • (6-3) 



1 — zze~ l ) \ (zz — e^t) (1 — zze^t) 

Proof. We start by inserting the expression ( |5.7| | of the operator e~ fLff into 
the equation ( |6.2| in which we also replace the GNBS by their definition in 
( |6.1) . We obtain successively 

Qm (L„) (z) = £exp(-fA£) <z,<r,m,ci \ tpt >< < | z,a,m,cc > (6.4) 

+00 

= ^exp(-U^) |<z/,ffl,« | < >| 2 (6.5) 



ir=0 



= £ exp (-fA£) (jV^ m (*))-* *«" (z) (6.6) 



fc=0 



= (A^, m (z))- 1 £ exp (-t (lot + k + 1)) (z) | 2 (6.7) 

fc=0 

= (A^, m {z))- l e-^ + ^ £ Or (z). (6.8) 

fc=0 

Now, to calculate the sum in ( |6.8| we make use of the expression of the 
functions 0^' m (z) in ( |3.8| l involving Jacobi polynomials and we apply the 
identity (4.8). This allows us to obtain the expression 

Qm (U) GO = (z))- 1 ,-^^ f ^"'" ^-f^ Y (6.9) 

V (1-zz) 2 / 

1-zz V ^ ^ 2^(1 -zz) 2 ^ 



l — zze-t) m I (zz — e~ f ) (1 — zze~ f ) y 

Finally we replace the factor A/"cr,m (z) by its expression in ( |4.10| . ■ 

In the following we will use the Q m — function presented above in order 
to write an inequality involving the thermodynamical potential associated 
with the operator L a . This potential reads 

O a := ^-Tr (Log (l + e"^-*)) ) (6 . 10 ) 

where r/ is the chemical potential and /3 = 1 / kg T, kg is the Boltzman con- 
stant and T denotes the temperature. Let us put e = e^ > and state the 
following inequality. 
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Proposition 6.3. Let a > 1. T/zen, f/ie thermodynamical potential in \6.10\ sat- 
isfy the inequality 



< n a (6.ii) 



max 

m6Z+n[0,(o--l)/2] 

/or epery /3 > 0. 

Proof. The form of the potential Cl a in ( 6.10| l suggests us to consider the 
function 

(p e (u) = -Log(l + eu). (6.12) 
So that we can rewrite ( |6.10| as 



tn K = Tr(<p 6 {e- a *)) t (6.13) 

where /3 = t G 1R+. We now apply the Berezin-Lieb inequality ( |2.7| | for the 
lower symbol Q m of the operator e~ fL « in the GNBSs representation ( |6.2| to 
obtain the following inequality 

J [cp £ o (L a )] (z) d ¥(r , m (z) < Tr {<p e {e- tL «)) . (6.14) 

Making use of ( |6.12| | and replacing the right hand side of ( |6.14| | by tO a as 
in ( 6.131 1, we get an inequality that holds for every m = 0,1, [(a — 1) /2] . 



Therefore, we consider the maximum with respect to the integer m of the 
quantity in the left hand side of ( 6.14| in order to be close as possible to the 



value of Cl x . 
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